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Abstract 

The problem of helicopter mechanical insta- 
bility is considered for the case where one blade 
damper is inoperative- It is shown that if the hub 
is considered to be nonisotropic the equations of 
motion have periodic coefficients which cann ot be 
eli mi nated- However, if the hub is isotropic the 
equations can be transformed to a rotating frame 
of reference and the periodic coefficients elimi- 
nated. The Floquet Transition Matrix method is 
shown to be an effective way of dealing with the 
nonisotropic hub and nonisotropic rotor situation. 
Time history calculations are examined and shown 
to be inferior to the Floquet technique for deter- 
mining system stability. A smearing technique 
used in the past for treating the one damper inop- 
erative case is examined and shown to yield uncon- 
servative results. It is shown that instabilities 
which occur when one blade damper is inoperative 
may consist of nearly pure blade motion or they 
may be similar to the classical mechanical 
instability. 
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The problem of mechanical instability of 
helicopters on the ground has been recognized and 
understood for many years. The analysis by Coleman 
and Feingold 1 has became the standard reference on 
this phenomenon although it was not published until 
many years after the first incidents of mechanical 
instability, or ground resonance as it is commonly 
known, were encountered on the early autogyros. 

The mechanical instability phenomenon is most com- 
monly associated with helicopters having articu- 
lated rotors; however, helicopters using the soft- 
inplane hingeless rotors which have become popular 
in recent years are also susceptible to this 
problem. Machines employing these soft- inplane 
hingeless rotors are also known to experience a 
similar problem, commonly known as air resonance, 
which occurs in flight rather than on the ground, 
fphfa air r^snnance problem has received much atten- 
tion in recent years (see, e.g.. Refs. 2 and 3) • 


will yield correct results, however, provided the 
equations are integrated over a sufficiently long 
time period. 

The purpose of this paper is to present a 
method of obtaining the mechanical stability char- 
acteristics directly for a helicopter operating on 
the ground with one blade damper inoperative. As 
will be shown later, the equations governing the 
motion of this system have periodic coefficients. 
This fact suggests the use of Floquet theory as the 
means for determining the stability characteristics 
of the system. In the following, the one-damper- 
inoperative problem is formulated and the resulting 
equations are solved using the Floquet Transition^ 
Matrix method described by Peters and Hohenemser. 
Results obtained using this method are compared 
with results obtained from the two previously used 
methods and recommendations are made concerning 
the future use of the three methods described. 


From the analysis of Reference 1 and others 
it is known that the ground resonance problem is 
due primarily to a coupling of the blade inplane 
motion with the rigid body degrees of freedom of 
the helicopter on its landing gear. These analyses 
have shown that with the proper selection of blade 
lag dampers and landing gear characteristics the 
problem of mechanical instability can be eliminated 
within the operating rotor speed range. All of the 
mechanical instability analyses conducted to date 
have one assumption in common - all blades are 
assumed to have identical properties. This is a 
reasonable assumption under ordinary circumstances; 
however, the U.S. Army has a requirement on new 
helicopters which invalidates this assumption. The 
requirement is that the helicopter be free from 
ground resonance with one blade damper inoperative. 
As will be shown later, this one blade damper inop- 
erative requirement has a serious impact on the 
classical method of analyzing a helicopter for 
mechanical instability. Further, there is at pres- 
ent no published method available for treating the 
case where each of the blades is permitted to have 
different properties. Thus the designer is faced 
with the dilemma of trying to satisfy the require- 
ment with an analysis method in which one of the 
basic assumptions is severely violated. 

Two methods have been used to circumvent this 
difficulty. The first of these involves a physical 
approximation so that the classical analysis 
becomes applicable. In this approach all blades 
are still assumed to have identical lag dampers 
even when one blade damper is removed, but the 
value of each of the dampers is reduced by the 
amount ci/N where N is the number of blades 
and c± is the original lag damper rate. As can 
be seen, with this approach a system is analyzed 
which is quite different from the actual situation 
of a rotor with no damping on one blade. The sec- 
ond method which has been used is to reformulate 
the equations of motion allowing for differing 
blade characteristics and to obtain the stability 
characteristics of uhe system using a time history 
integration of the equations. This second approach 
has the drawback that interpretation of stability 
characteristics from time history calculations is 
often difficult and open to question. The method 


Equations of Motion 

The equations of motion for the mechanical 
instability problem will be formulated using an 
Eulerian approach. It will be assumed, as is done 
in Reference 1, that the helicopter on its landing 
gear can be represented by effective parameters 
applied at the rotor hub. It will be further 
assumed that only inplane motions of the hub and 
blades are important in determining the ground 
resonance characteristics of the helicopter. Thus 
the degrees of freedom to be considered consist of 
two inplane hub degrees of freedom and a lead- lag 
degree of freedom for each blade in the rotor. The 
mathematical model to be used in the analysis is 
shown in Figure 1. Note that in the figure only a 
typical blade is shown. The analysis will be 
formulated for a rotor having N blades, and each 
blade is assumed to have a rotational spring and 
damper which act about the lag hinge. 


The blade equations are developed by summing 
moments about the lag hinge. The coordinates of 
the elemental mass dm in the fixed system are 

x i = x h + e cos ^ + P cosC^ + ) 

y± = y h + e sin ^ + p sinC^ + ) 


where 


♦ = ftt + 2*(i - 1)/N i - 1,2,..., N 


These expressions can be differentiated twice with 
respect to time to yield the accelerations exper- 
ienced by the differential mass 

= x^ - eft 2 cos - p(ft + t ± ) 2 cos 

- p( ± sinU ± + £ ± ) 

y ± = y fa - eft 2 sin ^ - p(ft + ^ i ) 2 sin(t i + 5^) 

+ p s’i C06 ^i + £i) (2) 


Using D , Alembert T s principle the summation of 
moments about the lag hinge can be written as 
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/p sin (t t + ^S^dm - / P cos (t ± + S ± )5r ± am 

- Vi - Vi - 0 i ~ 1 , 2 ,..., N ( 3 ) 

where the integrals are evaluated over the length 
of the blade. Introducing the expressions for 
x i and y^ and defining the following 


w 


the blade equations become 

Vi + eQ2s b sin 6in( +i + V 

- y h cosd^ + q)] + = o ( 5 ) 

i = 1,2,..., N 

If small displacements are now assumed the blade 
equations may be linearized to obtain 

*?i + Vi + (“Oi + n 2 v 2 oK = sin ^ 

- y h COS t ± ] i=l,2,...,H (6) 

where the following parameters have been introduced 

v o - eS A>' 

\ = V x b \ 

ni - W 

Under the assumptions stated earlier the hub 
equations of motion can be written directly as 

Vh + c x*h + Vh = P x 


Vh + Vh + Vh = P y 


( 8 ) 


where the coefficients on the left side of these 
equations are the effective hub properties in the 
x- and y-directions, respectively. The determina- 
tion of these properties depends on an extensive 
knowledge of the helicopter inertial character- 
istics and the stiffness, damping, and geometrical 
characteristics of the landing gear system. These 
properties may be determined either by ground shake 
tests of the helicopter, as suggested in Refer- 
ence 1, or by direct calculations. The right-hand 
side of the above equations are the forces acting 
on the hub due to the fact that the rotor is 
experiencing accelerations in the x- and y- 
directions. If the accelerations of the rotor 
center of mass are x and Jr c , respectively, the 
P x and Py are given by 


-a-Vc 

-"Vc 


( 9 ) 


The equations as written also indicate that in the 
absence of the rotor the hub degrees of freedom are 


uncoupled. This is an approximation, but it is 
an assumption made in Reference 1 and one generally 
used in helicopter mechanical stability analyses. 

If all blades in the rotor are assumed to have 
the same mass distribution, the coordinates for the 
total rotor center of mass may be written as 


1 N 
= + I 2- X i 

i=l c 




( 10 ) 


where x^ c and y-j^ are the coordinates of the 
individual blade center of ms.ss 7 measured with 
respect to the hub. If the center of mass of the 
ith blade is a radial distance from the lag 

hinge 


x ic = e cos + P c cos(i|r i + 

y± c = e sin + p c sind^ + q) 

Making the observation that, for H > 1 
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R 

COS 

k=il 


t k - £ sin y = 0 
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the rotor center of mass coordinates become 
x c - V - ( P C / N) £ h sin \ 


y c = y h + (p c/ N) 2^ cos +i 


( 12 ) 


These expressions may now be differentiated twice 
with respect to time and the forces P x an d P 
obtained as y 

N 


P x = _Ifa b\ + S b E [(l'i- fl2 Ci)sin 1^ + 20^ cos t 

n r „ 

P y = _Nm b y h " ^ E V ± )cos y i -2J^ i sin 


( 13 ) 


The hub equations of motion thus became 

(ffl x + + C A + Vh = 

N r- 2 . l 

Sb E^^i - a ‘ +i + 2^. cos 


( 14 ) 


(m + 

v y 


+ Vh + Vh = 


w .. 

Z (t ± - ft ^cos if ± - 2n£ ± sin 

/ 

The equations of motion for the system thus con- 
sist of (N + 2) coupled second-order differential 
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equations with the coupling terms having periodic 
coefficients. The periodic coefficients arise 
because the blade equations are written in a rotat- 
ing reference system whereas the hub equations are 
in a fixed system. As is shown in the Appendix, 
if all the blades have identical lag springs and 
lag dampers, the periodic coefficients may be 
eliminated through the use of multiblade coordi- 
nates. The effect of these coordinates is to 
transform the blade equations from the rotating to 
the fixed system of reference. The resulting con- 
stant coefficient system of equations is the set 
normally solved in the classical ground resonance 
analysis. As is shown, however, if the blades are 
allowed to have different lag springs and dampers, 
the periodic coefficients cannot be eliminated in 
the usual manner. 


y + + (a£ - Cl 2 ) y + 2 foe + Ot^x 

= E sinfiO-l) 

( 17 ) 

where the following parameters have been introduced 

v l = V (m x + Nm b ) ) 

= k x /( m x + 

= ° X / (n Sc + Nm b ) ) 

Introducing the rotating coordinates into the blade 
equations, Equations (6), results in 


An alternative does exist, however, for 
eliminating the periodic coefficients even when 
the blades are allowed to have differing character- 
istics. The alternative consists of transforming 
the hub equations into the rotating system of 
reference. In order to eliminate the periodic 
coefficients using this approach, the additional 
assumption must be made that the hub is isotropic. 
That is 
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This is the approach used in Reference 1 for treat- 
ing the two-bladed rotor which is another case 
where the periodic coefficients in the equations 
of motion cannot be eliminated by transforming the 
blade equations to the fixed system. 


The transformation from fixed to rotating 
coordinates is given by 

x « x^ cos at + y^ sin 0t j 

y = -x^ sin fit + y^ cos fit J 


( 15 ) 


Differentiating these expressions allows the fol- 
lowing identities to be established 

cos at + y^ sin Ot = x - fiy 

-x^ sin Ot + y h cos Ot = y + Ox 

1* 2- a 

x^ cos Ot + y^ sin Ot = x - fi x - 20y 

- 2 - - 

-x^ sin Ot + y h cos Ot =* y - fi y + 2 Ox 

The hub equations in the rotating system are then 
obtained by appropriate combinations of the x^ 
and y h equations, Equations (l^). The resulting 
equations are given below 

x + V' + - n )x - 20 y - Or^y 



h + + ( %j + 


= (vj/e) 


(x - a^x - 20y)sin (j - l) 


(y - fipy + 20x) cos ( j - l) 


at 


( 19 ) 


j = 1,2,..., N 

Since modern helicopters do not in general 
have isotropic hubs, the above equations can only 
be used to approximate the effects of a noniso- 
tropic rotor. They are, however, easily solved for 
the stability characteristics of the system and 
thus they might be used to obtain a first approxi- 
mation to the mechanical stability boundary for a 
helicopter with one blade damper inoperative . 


From the foregoing discussion it can be seen 
that if either the rotor or the hub is isotropic, 
the mechanical stability characteristics of the 
system may be obtained using conventional tech- 
niques. If both the rotor and hub are nonisotropic 
the equations of motion of the system contain 
periodic coefficients and thus the standard eigen- 
value techniques cannot be used to determine 
whether the system is stable or unstable . It is 
the purpose of this paper to demonstrate that 
Floquet theory can be used to analyze this general 
situation of a nonisotropic rotor coupled with a 
nonisotropic hub. 


Solution of the Equations 

If the periodic coefficients in the equations 
of motion are eliminated by assuming either an 
isotropic rotor or an isotropic hub, the stability 
of the system can be determined using standard 
eigenvalue techniques. The general case of a 
nonisotropic rotor coupled with a nonisotropic hub 
will be treated using Floquet techniques as 
described by Peters and Hohenemser , and Hohenemser 
and Yin. 5 A brief description of the technique 
will be presented here for the sake of completeness. 

In state vector rotation the free motions of 
the system may be written as 

J<5 j = [D(t)]|aj (20) 
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where the state variables for the problem being 
considered consist of 

• • • 

^1* C 2 > ^3J' y h' ^1' ^ 2 * ^K» *h* y h 

and the equations which describe the motions of the 
system are Equations (6) and (14). The matrix 
[D(t)] is periodic with period T and for the 
mechanical stability problem T = 2x/n. 

Floquet f s theorem states that the solution to 
the above system of equations has the form 

jaj -tA(t)]{ae^ +i ^ t | (21) 

where [A(t)j is the characteristic matrix and is 
also periodic with period T. The column of 
tial conditions jz(o)j is used in determining 
as 

|a j = [A(0)] -1 ja(0)| (22) 

The matrix [A(o)], the modal damping A, and the 
modal frequency cd are determined from the Floquet 
Transition Matrix [Q] which is defined by the 
equation 

jz(T)| = CQ]|a(o)| (25) 

for all sets of initial conditions jz(o)J . It is 
shown in References k and 5 that the eigenvalues 
Aj of the matrix [Q] can be used to determine Aj 
and o)j since 

Aj - (at) 

and the modal matrix of [Q] is just [A(o)J. The 
characteristic matrix [A(t)3 is then shown to be 
given by 

CA(t)] - [0(t)][A(O)]^- U+1<u) ^ (25) 

where the state transition matrix [0( t ) ] is defined 
by 

|a(t)( - [0(t)]|s(o)| (26) 

The characteristic multipliers Aj of the 
system are uniquely defined since the matrix [Q] 
is real; however, only the real parts of the 
characteristic exponents p. = + itOj are 

defined uniquely since * 

Pj » i ( In | A j| +1 arg Aj) (27) 

The imaginary part can only be determined within an 
integer multiple of 2*/T. This indeterminacy of the 
o>j causes no particular difficulty if one is only 
interested in the stability of the system. However, 
if one is interested in under standing the mechanism 
involved in any instability which might be found, 
this indeterminacy can be quite troubles cane . 

The Floquet Transition Matrix which is the 
basic element needed in the stability analysis is 
easily determined by a numerical integration of 
the equations of motion over one period T. If 
one desires to compute the characteristic functions 


[A(t)J the matrix [0(t)] is saved at each time 
point in the numerical integration to obtain [Q] . 
For the calculations of this paper, the fourth 
order Runge-Kutta method with Gill coefficients^ 
was used for the numerical integration. 

A comment is in order concerning the charac- 
teristic functions [A(t)]. The matrix [A(t)3 is 
a complex valued matrix and is determined at as 
many time points as desired. The computation of 
these functions can be relatively expensive and 
intepretation can be difficult. The interpreta- 
tion is made easier by the procedure outlined in 
Reference 5 far converting the complex functions 
into real functions which may be plotted as func- 
tions of time. The scheme used is essentially the 
same as that used when it is desired to plot as a 
function of time the modes of a system having con- 
stant coefficients. That is, for a conjugate pair 
of characteristic exponents 

p j - V ‘“j 

h ■ h - ‘“3 

the characteristic functions are also conjugate 
pairs. Thus the real modal function column for 
this conjugate pair of characteristic exponents 
will be given by 

jZj(t)} = J Aj(t) ) Aj+lfflj )t + 

(28) 

where jA,(t)j is the Jth column of [A(t)] and 
|Jj(t)j is the complex conjugate of this column. 

The purpose in performing these manipulations is 
to be able to plot the modal functions to deter- 
mine the relative magnitudes and phases of the 
various degrees of freedom in each mode. A discus- 
sion of this technique as it applies to constant 
coefficient systems is given by Meirovitch.? In 
this paper the exp(Ajt) is omitted from the above 
equation since it is simply a constant which multi- 
plies each component of the mode and causes each 
component to damp at the same rate. Thus the 
plots of the characteristic functions which are 
presented later in the paper will appear to be 
neutrally damped. 

In making the calculations for this paper it 
was found that the output from the calculation of 
the modal functions became so voluminous and these 
calculations became so expensive that the modal 
functions were only computed for selected points. 
Generally a sweep of rotor speed was made and the 
results examined. If an unstable region was indi- 
cated the rotor speed corresponding to the mayirmnn 
positive Aj was rerun and the modal functions 
calculated . 

Discussion of Results 

In order to demonstrate the application of 
the above-mentioned techniques and to obtain a 
general understanding of the effect of one blade 
damper inoperative on mechanical stability, a set 
of parameters were chosen. The parameters in the 
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mechanical stability analysis were chosen so as to 
be in the general range of interest for a single 
rotor helicopter and were such that the system was 
stable with all dampers functioning up to a rotor 
speed of 400 rpm. The parameter values chosen for 
the calculations are shown in Table 1. 

The parameters presented in Table 1 correspond 
to an isotropic rotor and a nonisotropic hub. In 
the following discussion results are presented for 
the case of an isotropic hub coupled with a non- 
isotropic rotor and a nonisotropic hub coupled 
with an isotropic rotor as well as the case of 
interest which involves a nonisotropic hub coupled 
with a nonisotrpie rotor. When an isotropic hub 
is mentioned, this means that the hub parameters 
in both the x- and y-directions were assigned the 
values shown in Table 1 for the x-direction. An 
isutrOpie ruuur implies that all dampers are 
operational and a nonisotropic rotor is meant to 
indicate that the lag damper has been removed from 
blade number 1. The analysis has been formulated 
in such a way that any number of blade lag dampers 
or lag springs may be removed to make the rotor 
nonisotropic. The results presented here, how- 
ever, only involve the removal of the lag damper 
from one blade. 

The case of an isotropic hub was first run in 
an effort to become familiar with the nonisotropic 
rotor results before proceeding with the more 
complicated Floquet analysis. The isotropic hub 
permits the equations to be transformed into the 
rotating reference frame and results in a system 
of equations with constant coefficients, Equa- 
tions (l6), (17), and (19), even with a noniso- 
tropic rotor. 

Figure 2 shows the results of the calculations 
for the isotropic hub with all blade dampers work- 
ing. Note that since the equations were solved in 
the rotating system, the frequencies in the lower 
portion of Figure 2 are plotted in the rotating 
system. The numbers attached to the different 
modes in Figure 2 and in subsequent similar figures 
have no significance other than to provide a label 
for the various modes. In Figure 2 the dashed 
lines represent the uncoupled hub modes . The 
uncoupled rotor inodes follow along the curves 
labeled 1,2 which also represent, in the terminol- 
ogy of Reference 5 , the rotor collective modes. 

Note that the uncoupled blade frequencies are zero 
for rotor speeds less than about 65 rpm. This is 
due to the fact that the blades are critically 
damped for these low rotor speeds. At the higher 
rotor speeds modes 3 and 4 are essentially rotor 
modes and modes 5 and 6 are essentially hub modes. 
At the lower speeds, however, due to the coupling 
between rotor and hub, mode 4 changes to a hub 
mode and mode 3 changes to a blade mode. Note 
from the damping plot that all the modes indicate 
stability over the entire rotor speed range. 

The results for one blade damper inoperative 
and an isotropic hub are plotted in Figure 3. 

Note that the removal of a blade damper has caused 
the appearance of a mode which was not present in 
Figure 2, namely the mode labeled 3 in Figure 3, 


and that this mode exhibits a mild instability 
between 160 and 200 rpm. At rotor speeds below 
about 100 rpm this mode has a frequency which cor- 
responds to the uncoupled frequency of the blade 
which has no damper. At rotor speeds above 100 rpm 
this mode begins to deviate in frequency from the 
uncoupled frequency. Another interesting point is 
that mode 1 in Figure 3 is precisely the same as 
the collective modes of Figure 2, and in Figure 3 
there is only one such mode. Thus it appears that 
the unstable mode in Figure 3 has evolved from one 
of the two collective modes shown in Figure 2 
because of the removal of one of the blade dampers. 

A time history calculation was made for the 
point of maximum instability in Figure 3 which 
occurs at approximately 175 rpm. The results of 
the time history calculation are shown in Figure 4. 
These results were obtained using the same inte- 
gration scheme as that used for generating the 
Floquet Transition Matrix. The top portion of the 
figure represents the individual blade lag motions 
whereas the lower portion represents the hub 
response in the x- and y-directions. Note from 
the figure that each of the degrees of freedom was 
given an initial displacement but the initial 
velocities were zero. The equations were inte- 
grated for 17 rotor revolutions. The figure indi- 
cates the blades which have lag dampers are well 
damped, but the blade on which the damper is 
inoperative experiences large lag excursions. 

Also, the hub motions, although not large, do not 
appear to have a high degree of damping. From the 
time history one would conclude that the system is 
stable since the motions of the various degrees of 
freedom do not appear to be increasing in ampli- 
tude with increasing time. The eigenvalue analysis 
has shown, however, that an instability exists. 

The problem with the time history calculations is, 
of course, that the equations of motion have not 
been integrated over a sufficiently long time 
period for the initial condi ticns chosen. Herein 
lies the difficulty with using the time history 
approach for calculating the stability character- 
istics of systems. One can never be sure if a 
sufficiently long integration period has been 
used, and the choice of initial conditions which 
will minimize the integration time required is a 
trial and error process. It has been observed on 
an analog computer that for the ground resonance 
problem the choice of initial conditions has a 
strong bearing on the conclusion inferred from the 
time history traces. The time history integration 
is also much more time consuming on the digital 
computer than the eigenvalue analysis . The time 
to generate Figure 4 which is for only one rotor 
speed was much greater than the time required to 
generate the eigenvalue results for all of Fig- 
ure 3 . It is thus concluded that whenever it is 
at all possible the eigenvalue approach to sta- 
bility calculation is to be desired over the time 
history approach. 

Having examined the case of one blade damper 
inoperative an an isotropic hub, the next logical 
step is to examine the more realistic situation of 
a nonisotropic hub. Before examining the one 
damper inoperative situation it was first desired 
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to confirm that the system was stable with all 
dampers working. The modal damping and frequency 
of the various modes with all dampers working and 
a nonisotropic huh are shown in Figure 5- As can 
he seen from the damping plot, all the modes are 
stable. In this case the equations of motion are 
solved in the fixed frame of reference and hence 
the frequencies are plotted in this frame. The 
dashed lines on the frequency plot represent the 
uncoupled system: the horizontal dashed lines 

being the huh modes and the slanted dashed lines 
being the rotor modes. Note that because the rotor 
modes become critically damped at low rotor speeds 
the two uncoupled rotor frequencies ccme together 
before reaching the origin. The uncoupled rotor 
lines also represent the collective inodes for the 
rotor. These modes are completely uncoupled from 
the other inodes and hence are not included in the 
eigenvalue analysis of the nonisotropic huh 
coupled with an isotropic rotor. The damping for 
the collective modes is exactly the same as that 
shown for modes 1,2 in Figure 2. 

The validity of the Floquet analysis was 
verified by comparing results from this analysis 
with results from both the rotating system analysis 
(isotropic hub) and from the fixed system analysis 
(isotropic rotor). In each case the results from 
the Floquet analysis were identical to results 
from the other analyses. 

Having thus established the validity of the 
Floquet analysis, results were obtained for the 
nonisotropic hub and one blade damper inoperative. 
These results are shown in Figure 6. Note that 
these results are very much similar to those shown 
in Figure 5 except that, as was the case with the 
isotropic huh and one blade damper inoperative, 
there are additional modes introduced. Also indi- 
cated is a relatively strong instability between 
210 and 305 rpm. The frequencies of the addi- 
tional modes which are introduced correspond, at 
low rotor speeds, to the frequencies of the 
uncoupled blade which has no damper. In the rotor 
speed range where the instability occurs, however, 
the frequency deviates from the uncoupled value as 
indicated by the mode labeled 3- In this range 
and at higher rotor speeds the mode labeled 5 is 
nearer the uncoupled blade frequency. It thus 
appears that far this case the instability is more 
a coupled rotor hub mode than a pure blade mode as 
was indicated for the isotropic hub. 

This conjecture is further strengthened by an 
examination of the modal functions. The modal 
functions for a rotor speed of 255 rpm, which is 
the point of maximum instability, are shown in 
Figure 7- The functions are plotted over a time 
period corresponding to one rotor revolution. Note 
from this figure that blade 1, the blade without a 
damper, has a significantly higher contribution to 
the mode than the other blades. Also from the 
plot of huh response it can be seen that the par- 
ticipation of the lateral hub degree of freedom, 
which has the higher of the uncoupled hub fre- 
quencies shown on Figure 6, is considerable. It 
is thus concluded from Figures 5 and 6 that the 


one damper inoperative situation can lead to a 
classical mechanical instability. 

Time history traces for this same condition 
are shown in Figure 8. These traces show the same 
general trends as observed in the case of the 
isotropic huh, that is, a large response of the 
blade having no damper and moderate responses from 
the other blades and the huh degrees of freedom. 
Again the time history traces are inconclusive 
regarding the stability of the system. 

One of the methods used in the past for 
treating the one blade damper inoperative case 
involves a smearing of the total blade damping. 

The reasoning for this approach is as follows. 

If the rotor has N blades then the total damping 
available in the rotor is Ncj[ where ci is the 
damping on one blade. If one damper is removed, 
the total damping becomes (N - l)c^ . Thus, using 
this approach, each blade in the rotor would be 
treated as if it had a lag damper equal to 
Ci (H - 1)/K. 

After an examination of the preceding one 
damper inoperative results it would be expected 
that this approach would lead to unc anservat i ve 
results. This is due to the fact that the insta- 
bilities encountered in the previous results 
involved large motions of the blade which had no 
damper. The smearing technique results in damp- 
ing, which is not greatly different fran the 
original value, being applied to each blade and 
thus the true situation is not adequately modeled. 

To illustrate this method, the nonisotropic 
hub case was analyzed using the smearing approach. 
The results from these calculations are shown in 
Figure 9- Note that although mode 3 becomes 
lightly damped the system remains stable through- 
out the rotor speed range considered. The fact 
that mode 3 approaches instability is attributable 
to the fact that this mode was not heavily damped 
in the original calculations. A run of the iso- 
tropic huh case, where all the modes were origi- 
nally well damped, indicated that the smearing 
technique resulted in well damped modes for one 
blade damper removed. The smearing technique is 
thus not recommended for treating the one blade 
damper inoperative situation since it leads to 
unconservative results. 

Since one way for eliminating the classical 
mechanical instability is to increase the blade 
damping, it was decided to attempt this approach 
on the instability indicated in Figure 6. The 
approach was to leave the damping identically zero 
on one blade and increase the damping on the 
remaining three blades. The results of this series 
of calculations are shown in Figure 10 where the 
region of instability is presented as a function 
of blade lag damping and rotor speed. As can be 
seen from the figure, increasing the blade danping 
on three of the blades has very little effect on 
the stability boundaries when one blade has zero 
damping. This result was somewhat expected since 
from the previous calculations it was observed 
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Appendix 

If the rotor is considered to be isotropic 
the periodic coefficients appearing in the equa- 
tions of motion can be eliminated through the use 
of multiblade coordinates similar to those 
described in Reference 1. These coordinates 
essentially transform the blade degrees of freedom 
into a fixed reference frame. The transformations 
are given by 

N 

ij = £ q sin if 
i=l 

Li = 2 L 003 L 

Differentiating these expressions leads to 
the establishment of the following identities 

L sin L = L - n Li 
£. L 008 L = Li + n L 
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i 5i Eln = h - n \ - «n 

q 008 3 * s ii " n2 ^n + 2 ^i 


It can be seen from these identities that the 
transformation is made by multiplying the blade 
equations. Equations (6), by either sin or 
cos and adding the equations. Crucial to this 
operation is the ability to remove the and 

03g from the summations. ‘Phis can only be done 
if all the blades have identical lag springs and 
lag dampers. If cue or more of the blades have 
differing characteristics, the and/or 

cannot be factored from the summation and hence 
the identities above cannot be applied. Thus, if 
one or more of the blades are permitted to have 
different lag springs or lag dampers, the periodic 
coefficients cannot be eliminated using the pro- 
cedure described in this Appendix. 


If Equations (6) are first multiplied by 
cos and summed and then multiplied by sin 
and summed, the following equations are obtained 
after introduction of the identities (A2) 


5 n + 


Vn - 
+ °Vl = 


2 d - V 2 ) - ^Itu * 2 


- y h E cos 2 +i 

i=l 1 


\ E Sin cos *, 
. i=l 


*1 + Vi “ [ fl2(l ' v o> - «Sj*i " 2ai 

- = (v^/e)t £ sin 2 


II 


i-1 


N ^ 1 

- E^ sin cos i|r i | 


(A3) 


Making the following observations that for N > 2 


N 

E sin + cos ♦, = 0 
i=l 1 1 


E cos 2 = E sin 2 t = H/2 
i=l 1 i=l 1 


the equations become . 

' 6 h + Vn -f^ 1 - v o> - %]*n + 2°Si I 

* °Vl = -(^e)y h ( ^ 

?i + - [« 2 d - v 2 ) - - aofin J 

- | 

These two equations describe the rotor mot ions in 

the fixed frame of reference . In terms of the 
variables described by Equations (Al) the hub 
equations. Equations (l4) , become 

(m x + 1^)3^ + = Bjj. 

(m y + m b )y h + Cy y h + y h = -aj u 

The stability of the rotor-hub system can now 
be determined using Equations (A4) and (A5) which 
have constant coefficients. This set of equations 
or a set similar to it is the one normally used in 
helicopter mechanical stability analyses. 

As a final observation, note tbat if. the 
blade equations. Equations (6), are simply summed, 
the following equation 

Vo + Vo + (“bi + = 0 < A6 > 

is obtained, where 

£ 0 = E (A7) 

i=l 1 



This equation represents the rotor collective mode 
and it may be observed that this equation is com- 
pletely decoupled from the hub degrees of freedom. 
Hence, the collective mode cannot influence the 
stability of the system and it is therefore not 
normally included in the mechanical stability 
analysis. 
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TABLE 1. PARAMETERS USED IN THE SAMPLE CALCULATIONS 


Number of blades 
Blade mass, 

Blade mass moment, 

Blade mass moment of Inertia, 1^ 
Lag hinge offset, e 
Lag spring, k ± 

Lag damper, 

Hub mass , 

Hub mass, m y 
Hub spring, k x 
Hub spring, k y 
Hub darnper, c x 
Hub damper, c^ 


if 

6.5 slugs (94.9 kg) 

65.0 slug- ft (289.I kg-m) 

800.0 slug-ft 2 (1084.7 kg-m 2 ) 

1.0 ft (0.5048 m) 

0.0 ft- lb/rad (0.0 m-N/rad) 

5000.0 ft- lb-sec/rad (4067.5 m-N-s/rad) 

550.0 slugs (6026.6 kg) 

225.0 slugs (5285.6 kg) 

85000.0 lb/ft (1240481. 8 N/m) 

85000.0 lb/ ft (1240481. 8 N/m) 

5500.0 lb-sec/ ft (51078.7 N-s/m) 

1750.0 lb-sec/ft (25539.3 N-s/m) 



ROTOR SPEED. Hz 


Figure 2. Modal damping and frequencies for iso- 
tropic hub, all blade dampers working. Fre- 
quencies plotted in the rotating system. 
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Figure 3- Modal damping and frequencies for iso- 
tropic hub, one blade damper inopei^ative . Fre- 
quencies plotted in the rotating system. 
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Figure 5. Modal damping and frequencies for non- 
isotropic hub, all blade dampers working. Fre- 
quencies plotted in the fixed system. 
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Figure k . Time history calculations far isotropic 
hub, one blade damper inoperative, H » 175 rpm. 
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Figure 6. Modal damping and frequencies for non- 
isotropic hub, one blade damper inoperative. 
Frequencies plotted in the fixed system. 
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Figure 7. Modal functions for nonisotropic hub, 
one blade damper inoperative, H = 255 rpm. 
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Figure 8. Time history calculations for noniso- 
tropic hub, one blade damper inoperative, 

H =* 255 rpm. 



Figure 9. Modal damping and frequencies obtained 
for nonisotropic hub, one blade damper inopera- 
tive, using the smearing technique. 
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Figure 10. Instability region as a function of 
blade lag damping for the nonisotropic hub and 
one blade damper inoperative. 
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